4avai  Research  Laboratory 

Washington,  DC  i0375-5000 


AD-A237  428 


1 

i 

i 

I 


I 


NRL  Memorandum  Report  6772 


Vlasov  Simulations  of  Very-Large- Amplitude  Wave 
Generation  in  the  Plasma  Wakefield  Accelerator 


Jonathan  Krall,  Glenn  Joyce  and  Eric  Esarey 

Beam  Physics  Branch 
Plasma  Physics  Division 


June  17,  1991 


DTIC 


91-03581 


Approved  for  public  release,  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


PorfTi  Approved 
0MB  No  0704  0189 


•’•jDiiC  'fOO''’''-.  of  s  '  »»■ -ur  oer  '^soO^ve  ■r'€  \jQinq  ti'Wp  tof  reviewing  tnilfuniv^r^  iearrn.nq  e«'St'ng  Odta  vOyfif% 

g<lT^e'’'''q  if’C  tr'e  aata  needed  a^o  corro'eting  jno  'e'-'ey*i'nq  tre  .onection  of  information  S<»rfa  comments  r^oaroing  this  Durden  estimate  or  an,  other  asoect  of  tn.s 

cOHect'On  0*  i-'rn'marion  rciud'ng  suggestions  ’or  r».auc»ng  m-s  ouraen  to  •‘Vashington  MeadOuarters  Se'.nes  Directorate  for  information  Ooerations  and  Reports,  W  IS  jeffersori 
Oa»‘S  H.qn^a*  Suite  '204  .-Hir^tcn  * -i  222Q2-alC2  and  tc  tn®  o»f,,  a o*  Management  ana  Budget  Rsperaior*  Reduction  Project  (0 /04-0 *88)  Washington.  CC  20503 


I.  AGENCY  USE  ONLY  (Leave  b'jnk) 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

1991  June  17 


3  REPORT  TYPE  AND  OATES  COVERED 


Interim 


Vlasov  Simulations  of  Very-Large-Amplitude  Wave 
Generation  in  the  Plasma  Wakefield  Accelerator 


5.  FUNDING  NUMBERS 


E-AI05-83ER-40017 


6.  AUTHOR(S) 

Jonathan  Krall,  Glenn  Joyce  and  Eric  Esarey 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 
Naval  Research  Laboratory 
Washington,  DC  20375-5000 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL  Memorandum 
Report  6772 


9  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 

DOE 

Washington,  DC 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


12a.  distribution,  availability  statement 


12b  distribution  code 


Approved  for  public  release;  distribution  unlimited. 


13,  abstract  (Mstimum  200  words) 

Simulations  of  the  plasma  wakefield  accelerator  are  carried  out  by  following 
the  time  evolution  of  the  plasma  distribution  function  in  one  dimension  via  the 
Vlasov-Maxwell  equations.  Simulation  results  are  compared  to  numerical  solutions 
of  the  nonlinear  relativistic  cold  plasma  equations  and  to  previous  theoretical 
estimations  of  trapping  and  thermal  effects  on  plasma  waves.  It  is  found  that 
highly  nonlinear  wakes  are  obtainable  in  the  vicinity  of  the  driving  beam,  where 
the  thtirmal  velocity  spread  of  the  plasma  is  reduced.  In  this  region,  wake  ampli¬ 
tudes  can  significantly  exceed  the  expectations  of  relativistic  warm  plasma  models 
and  agree  closely  with  cold  fluid  theory.  In  all  cases,  however,  particle  trapping 
and  thermalization  due  to  particle  scattering  from  the  large  amplitude  plasma  wave 
reduce  the  wake  to  below  the  nonrelat ivistic  wavebreaking  limit  after  the  initial 
accelerating  peak. 
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VLASOV  SIMULATIONS  OF  VERY-LARGE-AMPLITUDE  WAVE 
GENERATION  IN  THE  PLASMA  WAKEFIELD  ACCELERATOR 


I .  INTRODUCTION 

1-3 

Novel  plasma-based  acceleration  devices  are  being  actively 

researched  due  to  their  ability  to  support  acceleration  gradients  in  excess 

of  10  GeV/m,  which  is  greater  than  two  orders  of  magnitude  beyond  those 

obtainable  in  conventional  linear  accelerators-  The  plasma  Wakefield 

2 

accelerator  (PWFA)  is  one  such  device,  wherein  a  moderate-energy  electron 
beam  drives  a  plasma  wave  which,  in  turn,  accelerates  a  high-energy 
electron  bunch.  This  process  has  been  demonstrated  experimentally.^  Of 
interest  in  this  scheme  are  limits  on  the  obtainable  accelerating  gradients 
and  transformer  ratios  in  the  plasma  Wakefield. 

The  transformer  ratio  is  defined  as 

R  e  E^/E  (1) 

where  E_  is  the  peak  decelerating  field  experienced  by  the  driving  electron 
beam  and  E^  is  the  peak  accelerating  field  in  the  wake.  The  physical 
significance  of  the  transformer  ratio  is  that  the  energy  gained  by  the 
trailing  beam  in  a  single  acceleration  stage,  AV,  is  given  approximately  by 
AW  =  RWq,  where  Vq  is  the  energy  of  the  driving  beam.  Theoretical  results 
for  linear  Wakefields  have  suggested  that  driving  beam  pulses  that  are 
symmetric  in  the  axial  dimension  produce  transformer  ratios  that  are 
limited  to  R  <  2.^  Further  results  showed  that  higher  transformer  ratios 
may  be  obtained  by  using  a  nonsymmetric  beam  pulse^  or  by  operating  in  the 
nonlinear  regime.^  The  difficulties  asscciated  with  generating  a  shaped 
beam  pulse  add  to  the  attractiveness  of  the  nonlinear  PWFA  approach.  The 
viability  of  the  nuiilineax  FWIA,  itowevei ,  uia>  "oi;  iimiied  by  the  effects  of 
plasma  temperature  and  trapping  of  plasma  electrons  by  the  large  amplitude 
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wave.  For  the  present  study  we  will  consider,  via  simulation,  generation 

of  nonlinear  plasma  waves  by  a  symmetric  beam  pulse. 

The  properties  of  large-amplitude  waves  in  cold  plasmas  have  been 

2  7-12 

studied  by  a  number  of  researchers.  ’  For  nonrelativistic  plasma 

waves,  it  was  found  that  the  peak  electric  field  was  limited  by  the 

9 

nonrelativistic  wave-breaking  field: 


E  «=  mv  .  «  /e  <  E  .  E  mcoo  /e 
max  ph  p  wb  p 


(2) 


where  m  is  the  electron  mass,  v^^  is  the  phase  velocity  of  the  plasma  wave, 

e,  assumed  positive,  is  the  elementary  charge,  c  is  the  speed  of  light  in 

2  1/2 

vacuum,  =  (41111^6  /m)  ,  and  n^  is  the  equilibrium  plasma  density.  As 

14  -3 

V  .  approaches  c,  E  s  E  .  .  Note,  for  instance,  that  at  n  =  10  cm  , 
ph  max  wb  P 

13 

E^jj  =  1  GeV/m.  Coffey  ,  using  the  nonrelativistic  Vlasov-Poisson 
equations  and  a  "waterbag"  distribution  function  for  the  plasma  electrons, 
showed  that  this  amplitude  is  reduced  for  a  warm  plasma.  In  this  case, 


„  ,  /  1  8  1/4  _  1/2, 1/2 

E  =  (mv  . «  /e)(l  -  tV  -  tV  +  2y  ) 
max  '  ph  p  3  3  ^  / 


(3) 


where  y  •=  3T/mVpj^  and  T  is  the  plasma  thermal  energy. 

For  relativistic  plasma  waves,  E  >  E^j^  is  possible. 


plasma,  the  limit  is 


8 


For  a  cold 


E  =  (mc«  /e)'rF’(Y  u  -  1) 
max  '  p  '  ph 


1/2 


(4) 


2  2  -1/2 

where  Ypj^  *  (1  -  ^  '  This  suggests  that  as  v^^  approaches  c,  E  >> 

Eyj^  is  possible.  As  in  the  nonrelativistic  case,  this  limit  is  reduced  by 

«  A  ^  •  t 

thermal  Katsouleas  and  hori  peifoimed  a  calculation  similar  to 


2 


that  of  Ref.  13  for  the  relativistic  case.  Again  a  waterbag  distribution 
function  was  used.  They  found 

^  N  -1/4,1  1/2  1/4. ,1/2 

^max  =  U  )]  ,  (5) 

which  is  assumed  to  be  valid  as  long  as  it  gives  a  smaller  value  than  Eq. 
(4),  as  is  the  case  for  parameters  of  interest.  An  alternative  approach 
was  taken  by  Rosenzweig^^,  who  calculated  the  saturation  amplitude  of  the 
plasma  wave  in  two  ways.  Firstly,  from  an  energy  balance  argument 
involving  loading  due  to  the  trapped  portion  of  a  Gaussian  velocity 
distribution,  he  found: 


E  =  (mc«  /e)(mc^/4T)^^^  . 
max  P 


(6) 


Secondly,  from  a  warm  fluid  plasma  model,  he  found  a  similar  expression: 


2  1/4 

E  e  (mco)  /e)(4mc  /27T) 
max  P 


(7) 


In  Eqs.  (6)  and  (7)  it  was  assumed  that  v 


=  Vj^,  where  Vj^  is  the  velocity 

2  2-1/2 

of  the  driving  beam,  Yj^  >=  (1  -  )  >>  1  and  Yj^  »  1/p.  These 

approximate  expressions  are  interesting  in  that  they  are  (approximately) 


independent  of  Ypj^-  In  a  later  calculation,  employing  a  three-fluid  model 


for  the  thermal  plasma,  Rosenzweig  found 
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2  1/4 

E  =  (mc«  /eXYuinc  /27T)^  ,  (8) 

max  '  p  ' '  b 

which  gives  results  that  are  similar  to  those  of  Eq.  (5)  for  parameters  of 
interest  here.  This  is  not  surprising  in  that,  with  both  the  three-fluid 
and  the  waterbag  models,  there  are  no  plasma  electrons  with  p  >  Pj.j.jj.t 


3 


where  p  . ^  is  the  thermal  momentum  in  the  three-fluid  model  and  the 
cri  t 

surface  momentum  of  the  waterbag  in  the  waterbag  model.  In  either  case, 
the  onset  of  trapping  is  quite  sudden. 

In  this  paper  we  present  simulations  of  plasma  waves  driven  by  an  non- 
evclving  electT-on  beam  with  velocity  v^^.  These  are  carried  out  by 
following  the  time  evolution  of  the  plasma  distribution  function  in  one 
spatial  dimension  via  the  Vlasov-Maxwell  equations.  Direct  simulation  via 
the  Vlasov  equation  is  an  ideal  way  of  examining  thermal  and  trapping 
effects  because  the  artificially  high  temperatures  associated  with  typical 
particle  simulations  are  avoided. 

In  following  sections  we  first  describe  the  Vlasov-Maxwell  equations 

for  this  system.  Ve  include  a  discussion  of  the  cold  plasma  equations  to 

which  we  will  compare  our  results.  Ve  then  present  simulations  in  the 

linear  and  nonlinear  cases,  including  simulations  of  the  so-called 

nonlinear  PWFa^  in  which  a  plasma  wave  is  driven  by  a  beam  of  length 

»  1,  where  X  -  2iiv  ./«  ,  to  provide  R  »  1  and  E  /E  .  »  1. 

p  ph  p  +  vb 


II.  ID  NONLINEAR  VLASOV  FORMULATION 

To  simulate  plasma  wave  generation  we  use  the  Vlasov-Maxwell 
equations: 


ii  ^  F  li 

3t  ym  3z  ^  ®  z  3p^ 


0 


(9) 


and 


3E 

=  -4ne(n  .  n^^  -  n^)  ,  (10) 
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2  22  1/2 

where  f(2,p^,t)  is  the  distribution  function,  y  -  (1  +  ^  ^ 

the  driving  beam  density,  is  the  electric  field  and 

n  =  Jdp^  f  .  (11) 

For  this  system,  Maxwell's  equations  reduce  simply  to  Gauss's  law. 
Changing  variables  from  (z,p^,t)  to  (C=ct-2,p^,T=t),  we  have 


li 

at 


+ 


(c 


.it  _  af  o 
ym  ac,  z  ap^ 


(12) 


and 


aE 

=  4ne(n  +  ~ 

The  simulation  is  carried  out  by  solving  for  E^(C,t)  from  the  current 

value  of  f(CjP2tT)  via  Eqs.  (11)  and  (13)  and  subsequently  updating  f  via 

Eq.  (12).  Initially  f  is  uniform  in  z  and  Gaussian  in  p^,  varying  as 
2 

exp(-p^/2mTQ) ,  where  Tq  is  the  equilibrium  plasma  thermal  energy.  Ve 

represent  the  driving  beam  as  a  fixed  charge  shape  that  propagates  with 

velocity  v^^  but  that  otherwise  does  not  evolve.  The  details  of  the 

simulation  code  are  discussed  in  the  appendix. 

For  comparison  to  our  simulation  results,  we  consider  the  cold-fluid 

limit  of  the  Vlasov  equation,  which  we  discuss  here.  For  a  non-evolving 

electron  beam  with  v^^  =  c,  steady-state  may  be  assumed  and  derivatives  with 

respect  to  x  may  be  neglected.  It  is  convenient  to  work,  in  terms  of  the 

2 

normalized  scalar  potential,  ♦  =  e*/mc  ,  the  normalized  electron  velocity, 

t  =  V  /c,  and  the  normalized  electron  momentum,  u  =  yP  ,  where  y  = 
z  z  z 
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2  1/2 

(1+u^)  .  In  these  variables,  the  relativistic  Vlasov  equation  (12) 

becomes 


where  f  =  f(C,u  ). 

z 

Analysis  of  the  characteristics  of  the  above  equation  indicate  that 

there  exists  a  constant  of  the  motion,  w  =  y  -  u  -  ♦»  which  is  the 

z 

normalized  electron  energy  in  wakefield  frame.  Hence,  any  distribution 
which  is  a  function  of  this  constant  of  the  motion. 


f(w)  =  f(Y  -  -  ♦)  , 


(15) 


is  a  general  nonlinear  solution  of  the  relativistic  Vlasov  equation.  In 

terms  of  the  normalized  momentum  u^,  the  distribution  f(w)  is  related  to 

f(u  )  by  du  f(u  )  =  dwf (w)  |du  /dw I ,  where  |du  /dw|  *  [1  +  (w+4>)”^]/2. 

Ve  now  consider  the  cold  plasma  limit.  This  limit  is  valid  as  long  as 

the  thermal  velocity  is  small  compared  to  the  trapping  width  of  the  plasma 

wave,  v^,  «  |v  -  v  I,  where  v  is  the  bulk  longitudinal  motion  of  the 
th  '  p  z '  z 

electrons  in  the  plasma  wave.  This  inequality  holds  for  plasma  wave 
amplitudes  sufficiently  below  wavebreaking.  The  cold  plasma  electron 
distribution  function  is  given  by  f(w)  =  nQ5(w  -  1),  i.e., 

f(w)  -  "^^(y  -  u^  -  ♦  -  1)  ,  (16) 


where  n^  is  the  ambient  plasma  electron  density.  Analytically,  thermal 
effects  may  be  included  by  choosing  a  more  appropriate  distribution 
function  f(v),  such  as  a  vaterbag  or  a  Gaussian  distribution.^^ 
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Using  Eq.  (16),  various  moments  may  be  calculated,  i.e.,  the  electron 


fluid  density  n(C)  =  Jdu  f(u  )  and  the  electron  fluid  velocity  6  = 

z  z  z 

n  ^ Jdu^(u^/Y)f (u^) •  One  finds 


n  =  n^d  +  (!  +  <♦.)  2j/2 


(17) 


and 


.  [1  -  (1  +  ♦)^1/[1  +  (1  +  4)^1  . 


(18) 


Using  Poisson's  equation,  the  self-consistent  nonlinear  equation  describing 
♦(C)  is 


2 

^(2nj^/nQ  +  (1  +  ♦)~^ 


11  , 


(19) 


where  Equation  (19)  describes  the  generation  of  nonlinear 

Wakefields  in  a  cold  plasma  by  a  non-evolving  beam  ^^^(C)  with  c  c.  A 

similar  equation  may  be  derived  to  describe  the  generation  of  nonlinear 

12 

Wakefields  by  an  intense  laser  pulse. 

The  cold  plasma  equations,  Eqs.  (17'i-(19),  were  previously  derived  by 
Rosenzweig^ ’ using  cold  fluid  theory  auid  were  subsequently  used  to 
analyze  plasma  Wakefield  generation.  Numerical  solutions  to  Eqs.  (17)-(19) 
will  be  compared  to  the  simulations  discussed  below. 


III.  SIMULATIONS  IN  THE  LINEAR  REGIME 

Simulations  in  the  linear  regime  showed  excellent  agreement  with 

14  -3 

theory  for  T^  <  5  keV.  Here,  we  considered  n^  =  2  x  10  cm  ,  such  that 


7 


where  p  and  v  are  the  average  momentum  and  velocity  as  given  by  moments 

z  z 

f(z,p^,T).  In  the  linear  runs,  T  increased  over  the  several  plasma 
wavelengths  of  the  simulation  region  but  with  minimal  impact  on  the 
results . 
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IV.  SIMULATIONS  IN  THE  NONLINEAR  REGIME 


A  series  of  simulations  were  performed  with  the  beam  profile  given  in 
Eq.  (20)  with  n^^^/n^  =  0.1,  0.2  0.3,  0.4  and  0.5,  =  0.24  cm  e  = 

100  and  Tq  =  19  keV.  Peak  electric  fields  and  corresponding  theoretical 
values  are  plotted  versus  in  Fig.  1.  Simulation  parameters  were  as 

in  the  linear-regime  runs.  As  an  example,  the  beam  density,  electric  field 
and  perturbed  plasma  density,  n^^  «=  n  -  n^,  are  plotted  in  Fig.  2  for  the 

'^bO^'^O  '  ase.  Figure  3  shows  the  corresponding  results  from  the  cold 

plasma  equations.  Note  that  in  our  plotting  convention  >  0  is 

accelerating. 

In  the  most  highly  nonlinear  cases,  n^^^/n^  =  0.4  and  0.5,  several 

interesting  phenomena  were  observed: 

(1)  Excellent  agreement  with  the  cold  plasma  equations  was  observed  up 

to  the  first  accelerating  E^  peak.  In  this  region  the  temperature,  given 
by  Eq.  (22-23)  and  plotted  in  Fig.  4  for  the  =  0.5  case,  drops  by  an 

order  of  magnitude.  Results  did  not  vary  when  the  initial  plasma 

temperature  was  decreased  to  Tq  =  5  keV  (see  Fig.  4).  For  reasons  of 
numerical  expense,  it  is  difficult  to  drop  the  temperature  further  in  the 
code . 

(2)  Particle  trapping  was  observed  at  the  initial  accelerating  peak  in 

E^.  This  is  shown  in  Fig.  5,  where  the  plasma  electron  distribution  is 

plotted  versus  p^  and  C-  Here,  a  small  portion  of  the  electron 

distribution  was  accelerated  to  the  simulation  boundary  at  p  /me  ■= 

z  ^  rnsx 

30.7,  where  the  boundary  condition  dictates  f=0.  The  loss  of  plasma 

electrons  at  the  p  «  p  boundary  changed  the  results  by  <1X  relative 

z  z  f  nisix 

to  a  case  in  which  p  /me  =  60.0  was  used.  A  later  run  with  p  /me  = 

' 2 , max  2 , max 

5.12  shoved  no  change  in  either  of  the  first  two  accelerating  peaks  in  E_^. 
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This  suggests  that  the  trapped  particle  distribution  is  not  a  significant 
load  on  the  vakefield. 

(3)  Relativistic  lengthening  of  the  wake  was  observed.  This  was  not 
in  agreement  with  Ref.  15,  however,  wherein  X  =  4i[2T^c/£0p.  The  difficulty 
is  that  trapped  and  nearly-trapped  particles  increase  beyond  the 
expectations  of  the  cold  fluid  model.  Here,  y  is  the  average  value  of  y  as 
given  by  a  moment  of  f(C>P2»T)* 

(4)  The  Wakefield,  E^(C)i  decreased  after  the  first  peak  to  a  value 
less  than  the  nonrelativistic  wave-breaking  limit  of  Eq.  (2).  This  low 
value  of  E^  was  accompanied  by  an  increased  plasma  temperature.  Again, 
this  result  did  not  vary  with  Tq. 

(5)  Results  were  independent  of  y^^  for  y^^  >  10.  At  y^  =  5,  a 
significant  increase  both  in  trapping  and  plasma  heating  was  observed. 
Note  that  y^^  =  5  is  highly  unphysical,  given  our  non-evolving  beam. 

Both  the  close  agreement  with  the  cold  plasma  equations  in  the  region 
near  the  beam,  where  the  wakefield  is  larger  than  might  be  expected  for  a  T 
=  19  keV  plasma,  and  the  sudden  thermalization  behind  the  initial  E^  peak, 
where  the  wakefield  amplitude  is  somewhat  reduced,  are  unexpected  results. 

The  close  agreement  between  the  simulation  and  the  cold  plasma 
equations  in  the  region  near  the  beam  may  be  explained  by  a  simple 
argument.  In  the  region  near  the  beam,  plasma  electrons  are  accelerated  in 
bulk  in  the  direction  opposite  the  phase  velocity  of  the  wave.  This 
acceleration  reduces  the  longitudinal  thermal  motion  of  the  plasma 
electrons.  The  resulting  cold  plasma  distribution  closely  reproduces  the 
results  of  the  cold  plc.saa  equations.  From  Eqs.  (22-23), 

T  -  P  -  Jdp^  f6r(  (24) 

Y  y  Y 
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where  Ay  =  Y  -  Y  and  we  have  assumed  that  n(C)  is  constant  o’'er  the 

accelerating  region  and  that  y  »  1  and  Ay  <<  Y  for  regions  in  which  f  >  0. 

A  similar  relation  between  thermal  velocity  and  y  has  been  observed  in 

19 

conventional  induction  accelerators.  This  effect  will  be  evident  in  the 
highly  nonlinear  runs  presented  below. 

The  physics  behind  the  initial  peak  in  E^,  where  thermalization  of  the 
plasma  electrons  is  associated  with  a  lower  value  of  the  wake  amplitude,  is 
less  clear  than  the  situation  near  the  driving  beam.  Thermalization  is 
apparently  caused  by  scattering  of  the  plasma  electrons  from  the  large 
amplitude  plasma  wave  and  may  be  associated  with  the  fraction  of  plasma 
electrons  that  are  nearly-trapped,  remaining  in  the  region  of  peak  density 
for  times  of  the  order  of 

P 

Evidence  for  this  concept  is  provided  by  a  run  in  which  we  set  =  5. 
Lowering  v^j^  allowed  a  significant  increase  (relative  to  the  y  >  10  cases) 
in  the  density  of  trapped  and  nearly-trapped  electrons  and  an  increase  in 
the  thermal  energy  in  the  wake  by  approximately  a  factor  of  two.  Further 
evidence  is  provided  in  Fig.  4,  wherein  a  lower  value  of  Tq  increased 
(slightly)  the  thermal  energy  in  the  wake.  We  speculate  that  the  lower  Tq 
provided  a  more  coherent  nonlinear  oscillation  and  stronger  scattering  of 
the  electron  distribution.  The  temperature  diagnostic  is  difficult  to 
interpret,  however,  in  that  T  is  lowered  in  the  decelerating  phase  of  the 
plasma  wave  by  the  bulk  acceleration  mechanism  described  above.  In  the 
accelerating  phase  of  the  plasma  wave,  T  is  artificially  increased  by 
trapped  and  nearly-trapped  electrons. 

If  we  take  T  -  50  keV,  we  can  somewhat  explain  the  low  value  of  E^ 
behind  the  first  peak  by  the  theoretical  saturation  limits  given  in  Ref.  15 
and  quoted  as  Eqs.  (6)  and  (7)  above.  These  give  E  /E  ,  =  1.26  and  1.10 

1113  X  W  D 
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respectively.  Equation  (5)  gives  the  much  higher  value  of  =  2.22, 
suggesting  that  elements  of  the  model  of  Ref.  15,  particularly  the  trapping 
process,  represent  physics  that  comes  into  play  before  the  limit  of  Ref.  14 
is  reached.  It  is  interesting  to  note  that  our  results  are  independent  of 
Yjj  for  Yjj  >  10  whereas  Eqs.  (4),  (5)  and  (8)  vary  with  Yj^-  This  suggests 
that  expectations  based  on  Eqs.  (4),  (5)  and  (8)  may  not  be  meaningful  in 
this  instance  as  long  as  one  assumes  v  =  Vj^. 


V.  NONLINEAR  PLASMA  WAKEFIELD  ACCELERATOR 

For  a  beam  pulse  with  n^^  =  np/2  and  >  ^p»  theory  suggests  that 
transformer  ratios  R  >>  2  may  be  obtained.  We  consider  a  beam  pulse  of 
the  form: 


(c-v.  )T  <  C  <  (c-v.  )T  +  L. 


otherwise 


The  nonlinear  cold  plasma  equations,  (17)-(19),  have  been  solved  for  this 
case,^  Defining  x  e  (1  +  and  1.  *  ca^ft  -  x(E.)  within  the  beam  is 
given  implicitly  by 


.  1/2,  ..1/2  ,  ,,  ,.1/2  ^1/2, 

?.  =  x  (x  -  1)  +  lnl(x  -  1)  +  X  ]  . 


Defining  e  2iiLj^/Xp  and  x^  e  it  can  be  shown  that  the  maximum 

fields  within  and  behind  the  beam,  respectively,  are  given  by 


(1  -  1/Xf)‘ 
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mc« 


(28) 


E  =  - -  1)' 

+  e  r 


such  that  the  transformer  ratio  is 


1/2 

R  =  E  /E  =  xt  . 
+  -  r 


Vlasov-Maxvell  simulations  of  the  nonlinear  plasma  vakefield 

accelerator  were  performed  with  parameters  similar  to  those  of  the  linear 

14  -3 

and  nonlinear  cases  above:  n^  2  x  10  cm  ,  Tq  =  19  keV,  nj^  given  by  Eq. 
(25)  and  Yjj  ■=  100.  In  this  case,  however,  we  use  *  0.83,  1.63  and 

2.79  such  that  Eqs.  (26-29)  give  R  =  2,  3  and  4  respectively. 

Simulation  results  are  summarized  in  Table  I.  Cold  plasma  model 
results  are  also  given.  As  an  example,  the  beam  density,  electric  field 
and  perturbed  plasma  density,  n^  =  n  -  n^,  are  plotted  in  Pig.  6  for  the 

Ljj/Xp  *  2.79  case. 

As  before,  we  observe  a  reduction  in  the  thermal  velocity  spread  of 

the  plasma  in  the  region  near  the  beam,  excellent  agreement  with  the  cold 

plasma  equations  up  to  the  initial  accelerating  peak,  thermalization  of  the 

plasma  behind  the  initial  peak  in  E^,  and  a  drop  in  E^  to  E^  <  E^j^ 

thereafter.  In  the  L./X  =  2.79  case,  the  temperature  behind  the  initial 

b  p 


E^  peak  was  diagnosed  by  smoothly  setting  E^ 


0  beyond  this  point  and 


allowing  the  oscillating  plasma  to  "settle".  Results  indicate  50  keV  <  T  < 
150  keV. 

Agreement  with  the  cold  plasma  model  up  to  the  first  peak  is  made 
clear  when  Fig.  6  is  compared  to  Fig.  1  of  Ref.  7.  In  fact,  the  “ 
2.79  case  shows  a  peak  electric  field  which,  while  down  by  6%  from  the  cold 
plasma  result,  exceeds  the  expectations  of  Eqs.  (5-8)  for  a  T„  =  19  keV 
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plasma.  This  surprising  result  is  caused  by  the  reduction  in  the  thermal 
velocity  spread  in  the  plasma  in  accordance  with  Eq.  (24)  as  is  observed  in 
Fig.  7,  which  shows  plots  of  T  and  y  for  this  case. 

Our  observation  of  <  E^j^  following  the  initial  E^  peak  is 
interesting  in  light  of  Eqs.  (2-8).  Again  we  find  relativistic 
oscillations  of  the  hot  plasma  limited  in  a  way  most  closely  corresponding 
to  Eqs.  (6-7).  Taking  T  =  100  keV  after  the  first  peak,  Eqs.  (6)  and  (7) 
give  ^n,3x^^wb  “  1*06  and  0.93  respectively,  in  reasonable  agreement  with 
Figs.  6  and  7.  As  before,  Eqs.  (5)  and  (8)  give  higher  values  of  E 

IHctX 

VI.  CONCLUSIONS 

Simulations  of  the  plasma  wakefield  accelerator  were  carried  out  by 
following  the  time  evolution  of  the  plasma  distribution  function  in  one 
dimension  via  the  Vlasov-Maxwell  equations.  Results  were  in  surprisingly 
good  agreement  with  numerical  solutions  of  the  nonlinear  relativistic  cold 
plasma  equations  in  the  vicinity  of  the  driving  beam,  where  the  thermal 
velocity  spread  of  the  plasma  is  reduced  by  the  mechanism  of  Eq.  (24). 
This  reduction  in  the  plasma  thermal  velocity  allowed  the  generation  of 
wake  amplitudes  that  exceeded  of  the  predictions  of  relativistic  warm 
plasma  models  as  given  in  Eqs.  (5-8). 

Thermalization  of  the  plasma,  however,  apparently  due  to  particle 
scattering  from  the  large  amplitude  plasma  wave,  limited  the  wake  to  E_  < 
E  ,  behind  the  initial  accelerating  peak.  The  thermal  energy  in  this 

VD 

region  was  observed  to  increase  with  the  density  of  nearly-trapped 
particles  (those  that  remain  in  the  region  of  peak  density  for  times  of  the 
order  of  w”^)-  Wake  amplitude  limits  behind  the  initial  accelerating  peak 

Id 


were  in  reasonable  agreement  with  Ref.  15.  This  agreement  may  be 
serendipitous,  however,  given  the  complexity  of  the  observed  phenomena. 

These  results  were  insensitive  to  for  y^^  >  10  and  to  the 
equilibrium  plasma  temperature  over  the  range  3  keV  <  Tq  <  20  keV.  Results 
were  also  insensitive  to  the  presence  or  removal  of  the  trapped  portion  of 
the  distribution  function. 

Finally,  the  nonlinear  PWFA  concept  of  Ref.  7  was  found  to  be  viable, 
within  the  context  of  the  one-dimensional  simulation,  up  to  the  first 
accelerating  peak.  Here,  a  transformer  ratio  R  =  3.78  was  demonstrated  for 
a  case  in  which  the  theoretical  value  was  R  =  4. 
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APPENDIX:  VLASOV-MAXVELL  SIMULATIONS  IN  ONE  DIMENSION 


To  model  beam-plasma  and  laser-plasma  interactions,  we  have  developed 
numerical  solutions  of  the  relativistic  Vlasov-Maxvell  system  of  equations 
for  implementation  on  the  Connection  Machine  at  the  U.  S.  Naval  Research 
Laboratory.  Fir  further  numerical  discussions  of  the  relativistic  Vlasov 
equation  see  Ref.  20. 

The  code  is  formulated  in  (CtP  »P„»P  jt)  coordinates  where  C  =  ct  -  z 
and  T  =  t.  The  Vlasov  equation  in  this  case  has  only  one  spatial 
dimension: 


ii  +  (c  _  ^)M.  +  ^  r(E  +  ^  =  0 

3t  '  Yin^3C  3p^^'  z  r®c  ymc  '  ’ 


where  f(C,p  jT)  is  a  reduced  distribution  function  for  the  plasma,  p  and 

Z  X 

Py  are  momenta  which  are  determined  from  canonical  momentum  conservation,  y 
2  2  2  221/2 

=  +  (Pv  Pv  P^  ^  1  ®v  determined 

X  ^  z  X  y  z  X  y 

from  Maxwell's  equations. 

Here,  Maxwell's  equations  are  formulated  in  the  Lorentz  gauge.  In 
terms  of  the  potentials,  A^,  A^  and  a  =  A^  -  i. 


1  ^  o3^tt 
C  g^2  ^ 


An(J^ 


pc)  , 


(A2) 


4nJ  , 

X 


(A3) 


(A4) 
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and 


1 

c 


■3C3t 


=  4np  , 


(A5) 


where  J  and  p  are  the  current  and  charge  densities,  respectively.  The 
fields  are  given  by 


z  ~c  3t  c  3t  * 


(A6) 


E 


X 


(A7) 


.  3A 
c  3t 


(A8) 


B. 


3C 


(A9) 


and 


B 


y 


(AlO) 


The  Connection  Machine  handles  these  equations  well  if  they  are 

differenced  explicitly.  We  split  the  Vlasov  equation  into  the  following 

.  21 
two  equations: 


3f 

3t 


0 


(All) 


and 


17 


(A12) 


For  the  first  of  these  equations,  we  do  a  simple  upwind  differencing.  For 

22 

the  second,  we  use  a  flux-corrected  transport  (FCT)  algorithm. 
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L. /X 
b  p 

R 

R  [theory] 

E  /E  , 

+  vb 

[theory] 

0.83 

2.09 

2.09 

1.81 

1.83 

1.63 

2.90 

3.00 

2.71 

2.82 

2.79 

3.78 

4.00 

3.63 

3.86 

Table  1.  Simulation  results  for  the  nonlinear  PUFA  runs.  Theoretical 
results  are  from  numerical  solutions  to  Eq.  (17). 
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temperature  T  (dashed)  and  average  gamma  y  (solid)  for  the 
.79  case  show  a  reduction  in  T  as  expected  from  Eq.  (2A). 


